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ABSTRACT 

Constant load creep data were obtained for four plastics: polyethyl- 
ine, polypropylene, nylon, and polycarbonate. Viscoelastic responses of 
these materials proved to be non-linear in both time and stress. An 
empirical creep equation, based on a time dependent retardation rate 
constant, was developed to characterize the creep behavior and good cor- 
relation was obtained for various stress levels. By assuming the de- 
veloped relationship to be an equation of state, creep behavior of simple 
composites of plastics was predicted with good results when viscoelastic 
responses of the components were similar. The importance of stress-time 
history is demonstrated by using the creep equation to estimate mechanical 
behavior in other modes of testing. The concept of a retardation spectrum 


and its application to the materials studied is considered. 


TABLE OF CONTENTS 


section 
Ly Introduction 
Ze Linear Viscoelasticity 
Dus Linear Viscoelastic Models 
4. Non-Linear Viscoelastic Behavior 
Dic The Experiment 
6. Selection of an Equation to Characterize Creep Behavior 
re Data Analysis 
8. Conclusions 
BIBLIOGRAPHY 
APPENDICIES 
A. Tabulation of Results 
Bi Computer Program for Determining Creep Parameters 
C.. Computer Program for Predicting Creep Behavior of 
Composites 
D. Prediction of Stress Relaxation Behavior from the 
Empirical Creep Equation 
E. Prediction of Stress-Strain Curves using the 
Empirical Creep Equation 
ne The Box Distribution Function for Retardation 
Spectra 
G. Equipment 


Page 
1S 
IES, 
he 
28 
30 
33 
42 
58 


oo 


61 
67 


78 


87 


97 


104 


119 





Figure 


LO. 
11. 
2. 
ie 
14. 
iD. 


RGF, 


ty 


FS . 


ee 
20. 


Bal 


Za 


TABLE OF ILLUSTRATIONS 


Boltzmann Superposition Principle 

Hookean Solid 

Newtonian Liquid 

Maxwell Element 

Voigt Element 

Response of Viscoelastic Elements to Stress 
Generalized Voigt Model 

Generalized Maxwell Model 

Creep Behavior of a Viscoelastic Material 
Four-Element Model for Creep Behavior 

Creep Test Specimen 

Tensile Test Specimen 

Clamped Composite 

Comparison of Creep Curve with Simple Exponentials 
Retardation Spectra for Simple Exponentials 


Creep Behavior - Spike Distribution of Retardation 
Times 


Retardation Spectra - Monodisperse Amorphous 
Linear Polymer 


Spectra for Polymer of Mixed Molecular Weight 
Fractions 


Box Distribution of Retardation Times 
Method of Approximating Retardation Spectra 


Extrapolation of Creep Data to Determine Instantaneous 
Strain 


Correlation of Parameter Al with Stress 


Page 
18 
Dall 
21 
Za 
Zz. 
22 
23 
ZS 
24 
Ze 
31 
31 
31 
a5 
35 


36 


38 


38 


39 
oo 


43 


45 


Figure Page 
23, Stress-Strain Curves at Various Strain Rates 46 


eo Zero Time Strain Extrapolations from Stress- 46 
Strain Data 


Z5., Determination of Parameters Ay and A. 48 
26. Determination of Parameters C and Q 48 
27. Correlation of Parameter A, with Stress 30 
28. Correlation of Parameter A, with Stress St 
29. Strain vs. Ln Time without Viscous Flow a 
a0, Strain vs. Ln Time with Viscous Flow a 
Bl. Strain vs. Ln Time for Polycarbonate Creep Data 53 
B2. Predicted Relaxation Behavior of Nylon - Short Term 88 
B3% Predicted Relaxation Behavior of Nylon - Long Term 89 
34. Graphical Method of Predicting Creep Curves from 98 
Stress Strain Data 
35 Predicted Stress Strain Curve for Nylon 100 
36. Predicted Stress Strain Curve for Polycarbonate 101 
a7. Box Distribution of Rate Constants 105 
BO Wedge Distribution of Rate Constants 105 
Bon Delayed Strain Fraction vs. Ln Time for a 110 


Single Retardation Time 
40. Sigmoidal Curve for a Box Distribution 110 


41. Examples of Poor Characterizations of Delayed LO 
Elastic Response 


ae. Limiting Cases of the Box Distribution 114 


aE). Approximations of the Response Function for 114 
the Limiting Cases of the Box Distribution 


44. Delayed Elastic Response Functions Fit to Creep 117 
Data 


Figure Page 


or Monomodal Retardation Spectrum Approximated by 118 
Box Distributions 


46. A Polymodal Distribution Function 118 
47. Lever Arm Loading Assembly iL 
48. Specimen Grips and Clamps 121 
49. Eleven Unit Creep Testing Apparatus 22 
50. Tripolitis Extensometer Attached to Specimen E22 
Dk. Clip Gauge Extensometer P23 
Eye Clip Gauge Extensometer Mounted on Specimen 124 
Soe Bridge Circuitry for Extensometer 125 
54. Dimensioned Sketch for Clip Gauge Stress Analysis 126 
mie) Theoretical Stress Analysis of Clip Gauges 127 





rmn<e~ se «x 4 


QQ om @ 


oA 


e 
t 


Ce 
ae Cc 
t 


J 


Ea Co) 


NOTATION 
shear stress 
shear strain 
viscosity coefficient 
shear strain rate 
shear modulus 
tensile modulus 
time 
engineering strain 
strain rate 
engineering stress 
stress rate 
initial stress 
stress at time "Ui " 
creep compliance 
current time 
natural logarithm base 
glass transition temperature 
free energy of deformation 
zero time strain (instantaneous elastic strain) 
delayed elastic strain 
relaxation modulus 
non-recoverable flow strain 
deformation 
deformation rate 


polyethylene 


PP - 
PC - 


polypropylene 

polycarbonate 

nylon 

retardation time constant 

rate constant for delayed strain response 

rate constant distribution function 

retardation time constant distribution function 
delayed strain fraction 

creep parameter representing zero-time strain 

creep parameter representing steady-state creep rate 


creep parameter representing magnitude of delayed 
elastic strain component 


time dependent rate constant in empirical creep 
equation 


time exponent in empirical creep equation 
stress coefficients in generalized creep equation 


stress exponents in generalized creep equation 


10 


ACKNOWLEDGEMENTS 
The author wishes to express his appreciation for the guidance 
and encouragement given by Dr. G. F. Kinney in the pursuit of this 
project. The assistance of Mr. W. D. Roberts in equipment assembly 
and data taking proved invaluable. He is also grateful for the sug- 
gestions of Dr. W. Tolles which helped in the data analysis. Finally 
he wishes to thank his wife Phyllis for patiently following this 


investigation to its completion. 


Pe 





i. Introduction. 

Composite materials with high strength-to-weight ratios show great 
promise. [3] Composites may be looked upon as intimate mechanical combina- 
tions of two or more materials, but differ from multiphase materials such 
as metallic alloys in that they are not homogeneous on a macroscopic scale. 
Composite materials are not new. Concrete and filled resins, for example, 
have been successfully employed for years, and currently a bewildering 
number of combinations of many thousands of materials is available. Of 
necessity, however, the work with composites to date has been directed 
towards development of materials to meet specific technological require- 
ments, and very little research has been carried out in the area of 
characterization of the behavior of composites. [1, 2] 

Physical properties often constitute the starting point for the 
materials selection process, and so characterizing mechanical behavior 
of composites is an important first step in a fundamental research study. 
Of particular interest is viscoelastic behavior, frequently manifested by 
the component materials of the composite structure. Deviations from 
true elasticity are negligible for engineering purposes for materials 
such as metals and glasses, but in polymeric systems, which are extensively 
employed in composites, behavior is dominated by viscoelastic phenomena. 
{2] Many conventional composite materials have mechanical properties 
that exhibit appreciable time dependence through creep and relaxation pro- 
cesses. [4] The stress analysis becomes quite complicated, but it must 
be performed if a viscoelastic composite is to be utilized in an effici- 
ently designed structure giving satisfactory performance in service 
environments. [5] 


In this investigation the main purpose was to characterize the creep 
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behavior of individual viscoelastic materials by analytic equations, and 
then obtain suitable correlations so that the behavior of composites of 
these materials could be predicted... In addition, these analytic expressions 
were used to estimate performance of materials in modes of testing other 
than creep, such as stress relaxation and constant strain rate uniaxial 
tension tests. 

For the experiment, four different thermoplastics (polyethylene, poly- 
propylene, nylon, and polycarbonate) were chosen, and composites were 
formed as layered sandwich structures. A simple composite was selected 
for this basic study to eliminate as many variables as possible, since 
stress analysis of polymeric materials is inherently complex, even if 
linear viscoelasticity is assumed. In all experiments the specimens were 
subjected to uniaxial tension, and it has been shown [6] for this type of 
loading that changing the type and amount of component interfacial con- 
tact has negligible effect on composite behavior. This investigation was 
carried out in a laboratory with controlled temperature of 23 + 1° C. 
and relative humidity 50 + 2% in accordance with ASTM Standards. [7] 


Laboratory set-up and equipment is reported in Appendix G. 
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oe Linear Viscoelasticity. 

For the ideal elastic solid, stress is proportional to strain when 
a load is applied in accordance with Hooke's Law. It is also independent 
of the rate of strain. For example, at stresses below the yield strength, 


the behavior of most metals is given by Hooke's Law: 


(2.1) = Gad 


where 7 applied shear stress 
‘© = shear strain 
G = modulus of rigidity 
For the ideal viscous liquid, stress is proportional to the rate of 
strain, but independent of the strain itself, as stated by Newton's 


Viscosity Law. [8] At low shear stresses most liquids and gases obey 


Newton's Law: 


(2.2) id = 2¥ 

where P= applied shear stress 

= shear strain rate 
7} = Newtonian viscosity 
coefficient 

The laws and equations above represent idealizations. If sufficiently 
precise measurements are made it is found that no real solid is truly 
Hookean and that no real liquid is exactly Newtonian in mechanical re- 


sponse. [9] To represent such deviations an equation can be written in 


the form: 


(2.3) = et 


For the ideal elastic solid, M% equals zero and the equation becomes: 


(2.4) Y= ae 


ilps 


where the constant B is the reciprocal of the shear modulus, Vo . For 
a perfectly viscous liquid, the exponent M™ equals one, and the equation 
becomes: 

a ee 


or 


® 


(295) By ieee 


where the constant [) is the reciprocal of the Newtonian Viscosity, /n ; 
For a real material the exponent /“™ can vary between zero and one, and 

is a measure of the relative importance of the elastic and viscous contri- 
butions to the mechanical behavior. For values of the exponent ™ be- 
tween zero and one, equation (2.3) indicates a time behavior in which 
stress depends on both strain and strain rate. A material such as this 
which shows characteristics intermediate between those of an ideal solid 
or ideal liquid is called viscoelastic. 

In some situations neither the strain nor the strain rate is direct- 
ly proportional to stress. An example of this complicated relationship 
is the plastic deformation of metals. In the absence of this complexity, 
and when both strain and strain rate are linear with stress as indicated 
by equation (2.3) a material is said to display linear viscoelastic be- 
havior. [2] 

Linear viscoelasticity might be expected with materials whose struc- 
ture is not changing during the time it is subjected to an applied stress. 
[11] For such linear viscoelastic behavior, a law known as the Boltzmann 
Superposition Principle applies. This principle is illustrated in Figure 
1, which relates to a creep test where a material is subjected to a con- 
stant stress, and the creep strain changes with time. For linear visco- 


elastic behavior, this strain is directly proportional to stress: 


16 


(2.6) E(e— Care 


where ace) , the creep compliance, is a function representing the time 
effects of the viscoelastic response. For a stress O, applied at zero 


time, the strain as a function of time is: 


(2.7) e(t) p= On SCS 


For some higher stress 0, + O, applied at zero time, the strain would be 


increased: 


(2.8) a), = (Gr regs ©) 


Now if Oo is applied at time t=O, and stress O©, is added at a 


later time t = Te the strain becomes: 


(2.9) c(t} = Goemiee o -Tt-F,) 


This is illustrated in Figure 1 where linear superposition of the strains 


at their elapsed times is indicated. The stress history may be general- 


ized by: 
= 
(2.10) ae) — 2% Met) 


For a continuous stress history starting at zero time, the Boltzmann 


Superposition Principle can be expressed in integral form: 


e=t 
(2.11) e(t) = i J(t-2) o(e) de 


8=0 


es 
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where © is current time during the loading history, where o(e) is 
written for AO(e) » and t is constant in the integration. The creep 
strain in equation (2.11) is therefore represented by a convolution in- 
tegral where O(e) is the "forcing function," and 3(t-e) is the "re- 
sponse function." 


In a stress relaxation test, a tensile strain is rapidly applied and 


the stress changes with time. For linear viscoelastic behavior, 


(2.12) oCh) ea net) 


Here the viscoelastic time behavior is given by E-(t) , called the re- 
laxation modulus. Equations analogous to (2.7) through (2.10) can be 
written for successively applied strains, so that for a continuous strain 


history, 


e=t 
j E.(t-9) €(e) ae 


er 


(213) oO (t) 


It has been shown by Volterra that equations (2.11) and (2.13) are 
consequences of each other. [11, 12, 13] These equations are also a 
statement of Boltzmann's Superposition Principle and a definition of 
linear viscoelastic behavior. A thorough treatment of the mathematical 
theory and implications of linear viscoelasticity can be found in refer- 


ence [14]. 


19 


oe Linear Viscoelastic Models. 

Mechanical models are not necessary to describe viscoelastic phenom- 
ena, but they have been useful in describing a combination of Hookean and 
Newtonian characteristics of materials. The elastic mode of stress re- 
sponse may be represented by a spring (Figure 2), presumed to behave in 
an ideal Hookean manner, i.e., strain proportional to applied stress and 
independent of strain rate. The flow mechanism is indicated by a dashpot, 
(Figure 3) assumed to be filled with a Newtonian fluid so that its de- 
formation is proportional to the time and to applied stress, and inversely 
proportional to the viscosity of the fluid, N . The behavior of a 
material displaying flow superimposed upon elasticity can be represented 
by a combination of a dashpot and a spring. The two basic models are the 
Maxwell Element (Figure 4) and the Voigt Element (Figure 5). The response 
of these units to an applied stress is illustrated in Figure 6. In order 
to represent a real material an infinite number of these elements may be 
required. Generalized Voigt and Maxwell models are shown in Figures 7 and 
8. Poincaré [16] indicates that a phenomenon representable by any one 
model is also representable by an infinite number of other models. If 
one dashpot in the generalized Maxwell model has infinite viscosity 
(7| = C2 ) and if one dashpot in the generalized Voigt model has zero 
viscosity (7) = © ) then the models become equivalent in their representa- 
tion of viscoelastic behavior. 

The creep behavior of most materials can be considered as the super- 
position of three stress responses. (Figure 9) The instantaneous elastic 
strain can be represented by a Hookean spring, and the non-recoverable 
viscous flow by a Newtonian dashpot. The Voigt Element is particularly 
useful in describing a delayed elastic response. The viscous element 


coupled in parallel with the elastic element represents a damping resistance 


20 


Figure 2 
Hookean Solid 


Figure 4 


Maxwell Element 
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Figure 3 


Newtonian Liquid 


Figure 5 


Voigt Element 


(b) 





a 


Figure 6 


Response of Viscoelastic Models to Stress 


(a) Hookean Solid 
(b) Newtonian Liquid 
(c) Maxwell Element 
(d) Voigt Element 
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Figure 9 


Creep Behavior of a Viscoelastic Material 
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to 


Figure 10 


The Four-Element Model For Creep Behavior 
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to the establishment of elastic equilibrium. [10] A corresponding mech- 
anism in a polymeric solid might be the uncoiling of long polymer chains 
from their minimum energy unstressed positions, which involves diffusion 
of chain segments and takes considerable time. From the equations re- 


presenting the components of the Voigt element, 


deg 


(3.1) o Eé€4 + Ze 


t 


The solution to this governing differential equation gives the strain- 


time relation: 


C32) Cu 


£(i-e*") 


The creep behavior of linear viscoelastic materials illustrated in 
Figure 9 may then be described by the four-element model of Figure 10. 


The resulting equation becomes: 


Oo Oo -t/p 
(3.3) ——.— + e(i-£ + 2-t 

E>. E, No 
where “J = ye) is the retardation time, the time required for 
the delayed strain to reach a value | — Vo. , or 63.2% of its final 


value. 

The four-element model has been used successfully by Gearhart and 
Kennedy [17] in their study of cellulose acetate butyrate plastics, and 
by Dienes [18] to correlate data on creep of polyethylene at elevated 
temperatures. Marin and Pao [19] have extended the four parameter equa- 


tion to account for non-linear stress behavior frequently encountered in 
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the delayed strain and viscous flow components for a simple plastic. 


Their equation (in its simplest form) can be written: 
O my Pent mM 
(3.4) €e= 2+ c0™(i-¢%)+ Gomt 


A modified version of this equation was developed to characterize the 


creep behavior of the materials studied in this investigation. 


a 


4. Non-Linear Viscoelastic Behavior. 

In a previous study [6] involving polymeric materials in the same 
loading ranges used in this investigation, marked non-linear behavior was 
encountered in the elastic, delayed elastic, and flow regions of stress 
response. Such deviations from linear viscoelasticity are not uncommon. 
[20] In a study of creep in cellulose acetate, Haward [21] found that the 
delayed elastic strain component increased by a factor of six when the 
stress was increased 30%. This non-linear behavior occurred at a compara- 
tively low stress of 500 psi. Similar stress non-linearities have been 
observed for glass-like amorphous polymers [22], polymers in the transi- 
tion region just below He [2], and polymers in the rubbery state [23]. 

In all cases where deviations from linearity occur, an increase in stress 
produces a more-than-proportionate increase in strain. 

Alfrey [10] demonstrated that for large values of strain, non-linear- 
ity can be expected. From the theory of elasticity the work required to 
distort a solid to a given strain, called the "free energy of deformation," 


can be represented by a power series: 
2 = A 
ay Flv) = at™ + os? + ca¥"...... 


The shear modulus G has a magnitude given by: 


aie 
(4.2) G = d ¥% 





For small values of strain only the quadratic term is significant and 
the modulus is a constant, independent of strain = At larger strains, 
when higher order terms become important, the shear modulus becomes a func- 


tion of x , and stress non-linearities result. A similar analysis 
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predicts non-Newtonian flow of viscoelastic materials at higher stresses. 

For glassy polymers which have high moduli and can support large 
stresses, non-linear behavior has been observed even when the strains are 
relatively small. [2] This suggests that internal structural changes are 
responsible for the deviations. For example, it is reasonable to suppose 
that the flow rate will be markedly different for polymer chains which 
are straight than for a coiled chain. Experimental measurements in this 
field are sparse, but some theoretical aspects have been discussed by 
Schwarzlin terms of transition-state kinetics. [24] 

Nonlinearity of viscous and elastic elements of a model create 
systems of nonlinear differential equations which are difficult to solve 
by ordinary methods. In this investigation an equation similar to that 
of Marin and Pao (Equation 3.4) was used to fit creep data for various 
polymers. By making some simplifying assumptions in the equation, and 
with the aid of a digital computer, satisfactory predictions of composite 


behavior were obtained. 


Ze 


ore The Experiment. 

The viscoelastic materials selected for the investigation were: 
high density linear polyethylene, isotactic polypropylene, nylon 6/10, 
and polycarbonate. All were obtained from commercial sources in the form 
of extruded sheet. General mechanical properties of these plastics deter- 
mined by tension testing, along with approximate values of other physical 
constants from the literature [15, 25] are tabulated in Appendix A. 

Creep test specimens were cut and machined to the dimensions shown in 
Figure 11. The thickness of the specimens was that of the sheet materi- 

al; approximately 0.060 inches for polycarbonate and nylon, and 0.192 inches 
for polypropylene and polyethylene. Specimens used for relaxation tests 

and tensile tests were made with the dimensions shown in Figure 12. 

All specimens were cut with the long axis parallel to the extrusion direc- 
tron, 

Creep tests were carried out on constant load creep machines in a 
laboratory with the temperature maintained at 23 + 1° Celsius and relative 
humidity controlled at 50 + 2%. Construction features of the machines 
and the method of gripping the specimen are discussed in Appendix G. 
Series of creep tests were run at different stress levels, ranging be- 
tween 10 and 80 percent of the nominal yield strength. Extension of the 
test specimens over a two-inch gauge length was measured with a Tripoli- 
tis type extensometer which could be read to + 0.0001 inches. After load- 
ing the specimens, readings were taken at elapsed times of 6, 12, 30, 60, 
90, 120 seconds, and then in increasing intervals from 1 to 10 minutes 
during the first hour, then every 30 minutes during the first day. The 
test was continued for 100 to 150 hours, with an average of 50 data 
points taken during the run. Similar test procedures were followed for 
the composites, with the component materials clamped together outside the 
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Figure ll 


Creep Test Specimen 


Figure 12 


Tensile Test Specimen 





Figure 13 


Clamped Composite 
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gauge length to insure uniform elongation between the clamps. (See Figure 
13) 

Uniaxial tension tests were run at the same temperature and humidity 
on an Instron testing machine. Specimens were pulled at a constant rate 
of deformation, and Tripolitis extensometers were used to measure the 
extension over a two-inch gauge length. Preliminary experiments measur- 
ing the strain at timed intervals indicated that the strain rate was es- 
sentially constant throughout the test. 

Stress relaxation tests were performed on the Instron tester in the 
following manner. Samples were strained rapidly at ten inches per minute 
to a predetermined load level, where a load limit microswitch automati- 
cally stopped crosshead motion. Graphical output of load relaxation 
versus time at constant strain was obtained from the recorder installed 
on A testing machine. Initial straining was essentially instantaneous; 
less than one second was required to reach the desired stress. Dura- 
tion of the relaxation tests was approximately thirty hours. 

A portion of the creep data for polypropylene and polyethylene was 
obtained in a previous study. Because this work [6] was done in the 
same laboratory utilizing identical test procedures and specimens cut 
from the same stock, no distinction is made in the data tabulation in 


Appendix A. 
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ee Selection of an Equation to Characterize Creep Behavior. 

The four-element model (Figure 10) discussed in Section 3 describes 
the creep process in three stages; an instantaneous elastic strain, a 
delayed elastic response, and non-recoverable plastic flow. While there 
is considerable theoretical justification for this description of creep 
behavior [1l1, 27], this simple form of the model has two serious short- 
comings. First, it fails to account for the non-linear phenomena fre- 
quently encountered in tests made at stresses above 25% of the yield 
strength. Improved fits with actual data can be achieved by modifying 
the basic equation (Equation (3.3)) to account for these stress anomalies. 
Such a modification is shown in equation (3.4), and this method was uti- 
lized for data fitting in this investigation. 

t 

0 ( - a 4 or : t 


oO 
— ee + ae 
<4. 3) G. Ex E, | BC. 


(3.4) e= F = od tre toad et 

The second limitation imposed by the four-element model is more funda- 
mental. The presence of only one Voigt Element implies that the delayed 
elastic response of a polymer is a single molecular mechanism with one 
characteristic retardation time. As pointed out by Alfrey [15], the coils 
and linkages of a polymer chain have a multitude of sizes and spatial ar- 
rangements. One would expect local kinks to be straightened out more 
rapidly than long range coils, and that elastic equilibrium would be 
achieved sooner if the chain segment had a favorable orientation with \the 


stress axis. The delayed elastic response of a viscoelastic material 


would then consist of several mechanisms with different retardation times. 
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The "memory" effects observed in high polymers are frequently cited as 
proof of the existence of more than one retardation time. Failure of a 
single Voigt Element to represent accurately the delayed elastic response 
in the materials used in this study is discussed in Appendix F. 

Figure 14 illustrates schematically the shape of a creep curve for 
a viscoelastic substance compared to a set of simple exponential curves 
Grscthe form € = K ( =o eT) , and shows the failure of a single 
time constant concept to represent actual data. Figure 15 shows the re- 
tardation spectra corresponding to the simple exponentials. In this figure, 
the term b represents the rate constant (= Ur ) and y(b) the rate 
constant distribution function. The creep curve of Figure 14 can be ap- 
proximated very closely by a finite number of simple exponentials appropri- 
ately weighted. The model representing the delayed elastic response would 
then be the generalized Voigt Model of Figure 7. Such a representation 
consisting of exponentials @ ; @) » and S is given in Figure 16. 
A generalized equation for (MW discreet relaxation times, with the in- 


stantaneous strain (A,) and the flow strain (A, x t) included would be: 


> 


bat 
) 


(6.1) e= A, + Aq t * As(\- a, 


nM 


where Qj; = ¢ (bi) , the discontinuous distribution function for the ie 
exponential, and 7s OOL = 1. Although this "spike distribution" model 
was not used in this study to correlate creep data, G. F. Kinney, ina 
parallel investigation at NPGS, has used the model to fit the same data. 
In his work, rate constants were arbitrarily assigned, then appropriately 
weighted using a least squares curve fitting technique adapted to a digi- 
tal computer. Two to six retardation times were used, and excellent fits 


were obtained, even with the simple two-spike distribution. 
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Comparison Of Creep Curve With Simple Exponential Curves 
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Retardation Spectra Corresponding to 
the Simple Exponential Curves of Fig. 14 


35 


— —_— = 
— 


i weighted Sum of exponentials O,@S, Fig.\4 


actual creep curve 


time 


b= Vy; 


q(b) 





Figure 16 


Representation of Creep Behavior By Spike 
Distribution of Retardation Times 
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Another method of representing delayed elastic behavior is by a 
continuous distribution of retardation times, ¢ ( b ). This concept is 
attractive from a theoretical standpoint if one considers delayed elasti- 
city as consisting of many molecular mechanisms whose rate constants form 
a quasi-continuous spectrum. Alfrey [10] develops a molecular theory of 
delayed response and arrives at a theoretical distribution of retardation 
times for a monodisperse amorphous linear polymer. This distribution, 
whose upper and lower limits can be determined from the theory, is shown 
in Figure 17. (The notation here is left in the form given by Alfrey) 
For a polymer consisting of two different molecular weights, Alfrey has 
deduced the retardation spectrum shown in Figure 18. Tobolsky, Andrews, 
and Dunnel [28] used a box distribution (Figure 19), to correlate stress 
relaxation data. The limits of such a box distribution may ve obtained 
graphically as shown in Figure 19 from a data plot of delayed elastic 
strain versus logarithm of time. 

The calculation of the retardation spectrum for an actual material, 
where the time constants extend over many decades of time, can be an in- 
volved process. A simple method [15] of approximating the spectrum is 


(6.2) [.(T) = S(T 7 


Here,[ J(t) — t/n ] corresponds to the delayed strain response, 

ea ee ) » and ae is the distribution of retardation times which 
can be converted to the rate constant distribution q(b). Figure 20 
demonstrates this graphically by plotting the slope of the delayed strain 
vs. logtime curve versus logtime. [It is apparent that if Figure 20 does 
represent the response of the four-element model of Figure 10, the approxi- 
mation gives a curve with a peak at the most probable J instead of a single 
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Retardation Spectrum for a Monodisperse 


Amorphous Linear Polymer (Ref, [10] ) 
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Figure 18 


Spectrum for a Polymer Consisting of 
Two Different Molecular Weight Fractions 
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Box Distribution Of Retardation Times 
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Method Of Approximating Retardation Spectra 
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spike. For real materials, however, the approximation is better. Refer- 
ence [27] reviews more accurate methods for estimating the retardation 
spectra. 

Analysis of the box distribution function is given in Appendix F. 

It is shown that this distribution is unsuitable for characterizing the 
creep behavior of the materials tested, as are other monomodal distribu- 
tion functions. 

In the absence of a mathematically tractable function to represent 
the retardation spectrum of the materials tested, the continuous spectrum 
concept was discarded. Instead, a simple empirical equation was developed 
by assuming that the rate constant associated with the delayed elastic 
response is a function of time. One might expect that a polymer chain 
segment should exhibit a changing retardation time as the uncoiling pro- 
cess continues and the shape of the molecule changes. Thus, an overall 
characteristic time constant 7R might be considered, with a simple 


power dependency on time: 
mM 
(6.3) a at 


Substitution of this expression into the equation for the four element 


model (Equation (6.4)) leads to an equation of the form of equation (6.5) 


(6.4) e= A, + At + Aa (\— Vm ) 


_ 4 
(6.5) Ee=A, + Art + A(\-o% : 
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ie Analysis of the Data. 

Creep data were taken for all specimens in the form of elongation 
versus time. These data were punched on cards and a print-out obtained 
from an IBM 360 Digital Computer. The printout served as a check to pick 
up gross errors due to equipment malfunction during the testing period or 
errors made in transcription of the data. 

The first step in fitting the creep data was to obtain the time- 
independent part of the elastic strain. This instantaneous elastic 
strain, Eo » occurring when the load is applied, is represented by the 
parameter ey in equation (7.1): 


CQ 
Cab E= A, + Bot * ss hilaeiatieee ) 


To evaluate the constant ay » a direct method was employed. The first 
four data points were used to perform an extrapolation back to zero time. 
This extrapolation to determine the parameter os is illustrated in 
Figure 21. It is clear from this figure that this evaluation of A 
is at best a rough estimate due to the nonlinear nature of the extrapola- 
tion. It will be shown in Appendix D that the parameter A, does in- 
clude a portion of the delayed elastic strain, but for the time frame of 
this experiment, the strain Ay occurs essentially instantaneously. 

To account for the non-linear viscoelastic behavior encountered, the 
parameters Ay ; A> » and A3 were correlated with stress so that 
the creep behavior could be predicted for any stress level. The general- 


ized form of equation (7.1) becomes: 


ts 


) 


ae C= Ko + Bort + Po“ (\-2~ 
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4 
Time , nours x\oO3 


Extrapolation of Creep Data to Determine Instantaneous Strain 


The magnitude of parameters Aa and Aa, depend on the value chosen for 
Ay . Prior to evaluating Bue and Ax3 » the extrapolated values 
of zero time strain cs were correlated with the first term of equation 


(7.2)? 
N 
(7.3) AY, = KO 


7.4) JmA, = InK + NWO 


Values of Im Ay for different stress levels were read into a digital 
computer, and equation (7.4) was used to perform a least squares fit to 
determine the stress coefficient K and the stress exponent N 
Figure 22 shows the plot of equation (7.4) and the fitted curve. The 
data displayed in this figure are that obtained from the extrapolations 
of Figure 21, made for polycarbonate. 

An independent method of obtaining the zero-time strain was attempted 
by running a series of tensile tests at various strain rates, and plotting 
the resulting data as sketched in Figure 23. From these stress-strain 
curves, a cross-plot of strain versus reciprocal strain rate may be con- 
structed as illustrated in Figure 24. The zero-time strain can then be 
obtained by extrapolating the data to strain rate E = OO » which 
corresponds to instantaneous loading. This technique was attempted using 
a digital computer to perform a parabolic extrapolation, but the exponential 
character of the data plot invariably resulted in high estimates of the 
zero-time strain, Eo , as indicated in Figure 24. In addition, the 
maximum deformation rate available on the Instron tester used for the 


tensile tests was % = 10 inches/minute. Consequently, there was a lack 


of data in the region of greatest interest, i.e., where => 0. This 
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Correlation of Parameter Ay with Stress 
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method may be expected to produce excellent results if high-speed testing 
apparatus were available [29], but in view of the limitations discussed 
above it was discarded in favor of the creep data extrapolation technique. 


After the parameter A, was found by extrapolation and "smoothed" by 


L 
the stress correlation (Equation (7.3)) the creep data, along with the 
calculated parameter Ay» was placed in the digital Sen pumer vad program 
FITDATA (See Appendix B). The viscous flow parameter A, was determined 
by a linear extrapolation using data points at elapsed times greater than 
70 hours. It was assumed that elastic equilibrium was achieved by this 


time. With parameter A, Specified and parameter A, computed, the maximum 


Z 
magnitude of the delayed elastic strain, A; is determined. This pro- 
cedure used by program FITDATA is outlined in the sketch of Figure 25. 
The time constant parameters were computed by FITDATA by evaluating the 
delayed strain fraction, €& , at various times. 


_ ct 
(7.5) a = | ae 


from which 


~ bw (| - €) = ct? 
(7.6) In |-dn(\- ©) = Inc + aula 


A plot of the left-hand side of equation (7.6) versus |) ee was given 
a least squares curve fit, and the constants C and Q determined. The 
method is sketched in Figure 26, and graphical computer output of actual 
data is displayed in Appendix B. 

The creep parameters A, and A. were correlated with stress using 


2 3 
equations (7.7) and (7.8): 


(7:87) n= mot 
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Determination of Parameters C and Q 


48 


_ 


M 
(78) Az = Pd 


The method used for the A, correlation was employed to determine the stress 


1 
exponents and coefficients in equations (7.7) and (7.8). These determina- 
tions for polycarbonate are illustrated in Figures 27 and 28. 

A question arises as to the validity of inclusion of the viscous 
flow term Art in equation (7.1). If the term is not included and 
the first and second derivatives of strain with respect to log time are 
taken, at long times the second derivative is always negative and becomes 
vanishingly small. At long times the first derivative, corresponding to 
the slope of the € vs. bnt curve of Figure 29, approaches zero. With 
no viscous flow, elastic equilibrium is eventually reached and extension 
ceases. If the flow term At i included, the second derivative 
changes sign and becomes increasingly positive at long times. This behav- 
ior is observed as an inflection and an upswing in the plot of € vs. tnt ; 
Figure 30. If viscous flow is a factor in creep behavior, therefore, a 
data Riot of € vs. nt should show the inflection and upturn indi- 
cated by Figure 30. In all the materials investigated this behavior due 
to viscous flow was observed, as illustrated by the example of Figure 31, 
justifying inclusion of the tern A\2* in equation (7.1). 

Values of the creep parameters determined for the materials tested 
are summarized in Table 2 of Appendix A. The constants C and Q are ap- 
pera: ly, independent of stress and are of the same magnitude for all 
materials tested. Standard deviations in strain were acceptable for all 
runs. The errors in computed strain were less than 5% in the region where 
delayed elastic behavior is significant, except at elapsed times under 
two minutes. Sample computer output from FITDATA listing point by point 
comparison of observed and computed strain is given in Appendix B. For 


typical values of the constants Cw=0O.8, and Q = 0.35 , the delayed 
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Correlation of Parameter Ao with Stress 
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Correlation of Parameter Ag with Stress 
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Strain versus Ln Time with Viscous Flow 
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strain is only 97% exhausted at an elapsed time of /0 hours. For the 


extrapolation to determine parameter A elastic equilibrium is assumed. 


9? 
Consequently, in the resulting "forced fit," after the other parameters 
are determined, the computed strain is approximately one to two percent 
low at long times. 

From the tests at higher stresses, it is apparent that elastic equi- 
librium is not achieved during the testing period. This produces a high 
estimate of the parameter A,» which in turn results in a low estimate for 
parameter A... At low stresses the flow rates are very small and the er- 
ror introduced is negligible. However, power dependence of A, on stress 
creates noticeable deviations at high stress levels. This deviation is 
clear in Figures 27 and 28 at stress levels of O = 4500 psi and O= 
5000 psi. In performing the stress correlations, data at lower stresses 
was weighted more heavily, which corrects the erroneous values of A, and 


Z 


A. . 


(See figures 27 and 28). 

The fitted empirical equation (7.2) provides excellent predictions 
of creep behavior on individual materials loaded at constant stress. 

When a simple composite sandwich of two materials is tested in creep how- 
ever, the stress in each component changes with time due to the differ- 
ent viscoelastic responses of each material. If the time responses were 
linear in stress, accurate theoretical predictions of composite behavior 
could be made. Reference [15] outlines procedures for calculating stress 
as a function of strain and time when the creep compliance J(t) is 
linear in stress. The Boltzmann Superposition Principle (Equation (2.11)) 
could then be applied to determine the creep behavior of a composite. 
From Table 2 in Appendix A, the stress exponents are greater than one, so 
the viscoelastic behavior of the materials tested is non-linear and the 


creep compliance 3 (oe is not independent of stress. Resulting stress- 


strain-time relationships become extremely complex. If a non-linear 


54 


viscoelastic "equation of state"’ could be developed, it might be too compli- 
cated for practical use, even with digital computer techniques. 

In order to predict the behavior of composites studied in this investi- 
gation, a simplifying assumption was made. The strain response of a 
material was assumed to be independent of its stress-time history. Equa- 
tion (7.2) then functions as an equation of state; if a stress and time 
are specified, then the strain is determined. From the physical restric- 


tions placed on the composite tested in creep, 


(7.10) Ea = Et = € 


where P represents applied load and the subscripts a and b are used to 
identify quantities pertaining to component A or component B of the com- 
posite. Letting a and b denote the relative fraction of component A 


and component B respectively, 


(7.11) Quechua en) 
(i. 12) Ce Oggerela0 eX aT 


Unsubscripted symbols above pertain to the composite structure as a whole. 


From equation (7.2), 


7 +&b 
2.13). ya kyoetaBarktt +R (Vs .) 


Qa 
(7.14) EG, = Ka ae + B.On*t + PCa, = i ) 


ao 


Equation (7.12) may be solved for Oa in terms of known composite stress 
Oo » Known volume fractions a and b, and the unknown stress One 
With this result substituted for Oa in equation (7.14), equations (7.13) 
and (7.14) may be equated according to equation (7.10). For a specified 
time, the corresponding value of Oy may be found. A series of values 
of OF and time can then be substituted in equation (7.13) to determine 
strain €, according to the equation of state assumption. By equation 
(7.10), strain €4 is the same as the strain in the composite. 

A digital computer program, program COMCUR, was used to predict compos- 
ite behavior in this manner. Sample output from COMCUR consisting of point- 
by-point comparisons of predicted and observed strain is found in Appendix 
C, along with graphical displays of this output. A summary of results of 
this prediction method are tabulated in Table 4 of Appendix A. When the 
viscoelastic responses are similar, the load distribution in the components 
of the composite does not undergo large changes during the test and pre- 
dictions are excellent. Such is the case for the polyethylene/polypro- 
pylene (PE/PP) composites. It should be noted that if the stress in the 
components were constant, the method would be theoretically correct, since 
equations (7.13) and (7.14) are valid for a constant stress test. Large 
errors in predictions occur with component materials of nylon and poly- 
carbonate. As indicated by the magnitude of the creep parameters given 
in Table 2, Appendix A, the viscoelastic response of these materials is 
quite different. Graphical output from COMCUR shows that the load dis- 
tribution for the NY/PC composites changes rapidly during the first few 
minutes of the test. 

The empirical creep equation (Equation 7.2) was used to predict stress- 
strain behavior and stress relaxation behavior for single materials. This 


analysis was done with a digital computer using the equation of state 
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assumption. An explanation of program RELAX, used to predict stress re- 
laxation behavior is given in Appendix D. Output from program SIGEPS, an 
estimate of stress-strain curves, is given in Appendix E. [In both cases, 


computed responses were considerably different from data obtained by experi- 


mene. 
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Sz Conclusions. 

The results of computer programs RELAX and SIGEPS demonstrate the 
complex nature of viscoelastic behavior. In tensile tests and relaxation 
experiments the stress changes with time and the constant stress equation 
gives poor predictions of responses. While the equation gives excellent 
estimates of creep behavior, the stress-time history cannot be ignored 
when predicting viscoelastic behavior in other modes of testing. 

The combined mathematical and experimental approach to understand- 
ing composite behavior is a powerful tool. An analysis of simple compos- 
ite structures provides guidelines for predicting behavior of more compli- 
cated systems where factors such as component distribution and interfacial 
shear stresses become important. The non-linear behavior of many visco- 
elastic materials, even at moderate stresses, leads to involved mathemati- 
cal relationships. A digital computer is necessary for analysis of these 
complex systems. 

Some opportunities for further investigation are: 

(1) Use of X-ray diffraction techniques to follow structural 
changes during testing in order to determine the relation 
of internal structure change to observed non-linear visco- 
elastic effects. 

(2) Fabrication and testing of more realistic composite materi- 
als to determine if the empirical equation developed here is 
general enough to characterize the behavior of these materi- 
als in creep. 

(3) Use of other modes of testing, such as stress relaxation and 
dynamic mechanical tests to assist in development of non- 


linear viscoelastic theory. 
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APPENDIX A 


TABULATION OF RESULTS 


Table 1 


Properties of Tested Polymers 





*Tensile *Tensile *Percent!] Melt. 


Strength Modulus 


Material psi psi x 10° 
POLY ETHY LENE 4000 220 380 
POLY PROPY LENE 4400 190 210 
POLYCARBONATE 8800 300 60 
| NYLON 5900 120 185 
oe 


*At strain rate = 0.012 in/in/min 
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APPENDIX B 
COMPUTER PROGRAM FOR DETERMINING CREEP PARAMETERS 
A Fortran IV computer program, compatible with the IBM 360 system, 


was used to determine the creep parameters A C, and Q in the 


ci 

equation: ota 
Ee= A, tArt +As(\-Z2™ ) 

The program, called FITDATA, is entered with the predetermined parameter 

Ay and the deformation-time data obtained from the creep tests. Computer 

output consists of the calculated parameters, tabulated values of observed 

and computed strain, and intermediate data used in the parametric determina- 

tion. An off-line plotting routine, called SUBROUTINE DRAW, is used to 

display the output (experimental data and fitted curve) in graphical form. 

This subroutine was also used to plot the data for determining parameters 

C and Q in the equation above. Least squares fitting was done by the 

library subroutine, SUBROUTINE LSQPOL available at the NPGS Computer 

Facility. Program FITDATA, along with typical tabulated and graphical 


output, is given on the following pages. The program, discussed briefly 


in Section 7, is annotated with comment cards. 
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APPENDIX C 
COMPUTER PROGRAM FOR PREDICTING CREEP BEHAVIOR OF COMPOSITES 
As discussed in Section 7, the prediction of composite behavior was 

made by assuming the analytic expression for creep of the components to 
be an equation of state. Input data to the Fortran IV computer program 
called COMCUR consists of the experimental creep data for the composite, 
the volume fractions a and b for each component, and the stress coeffi- 
cients and exponents for each component. These were computed from stress 
correlations of the creep parameters for individual components obtained 
from program FITDATA, For a series of times ranging over the testing 
period, the stress in each component is computed. This is done by a 
Newton-Raphson iteration since the unknown cannot be expressed explic- 
itly as afunction of the other variables. With the stress distribution 
determined, the composite creep strain, equal to the strain in each com- 
ponent, is computed using the equation of state assumption. Subroutine 
DRAW is used to display the following output: 

(1) Calculated creep curve of component A at the given stress 

(2) Calculated creep curve of component B at the given stress 

(3) Predicted creep curve of composite AB at the given stress 

for the specified volume fractions of components A and B 

(4) Experimentally observed creep data for composite AB 
The load distribution in the composite during the first 200 seconds of 
the test is also presented in graphical form. Tabulated output and graphs 
for a typical run are given on the following pages. The program COMCUR 


and the related function subprograms are annotated with comment cards. 
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COMPARISON OF THE CALCULATEO ANO EXPERIMENTAL VALUES OF CREEP STRAIN 

FOR COMPOSITE RUN AAl1 1S GIVEN BELOW: : 

TIME »HOURS EXPTL. STRAIN CALC. STRAIN PCTs ERROR 
0.00056 0.005500 0.004749 13.66 
0.00139 0.006200 0.005013 19-15 
0.00278 0.006400 0.005274 17.60 
0.00417 0.006650 0.005456 17.96 
0.00556 0.006800 0.005599 17.66 
0.00694 0.006900 0.005719 17-12 
0.00833 0.907000 0.005823 16-81 
0.01111 0.007150 0.006001 16.08 
0.01388 0.007300 0.006148 15.78 
0.01667 0.007400 0.006277 15.18 
6.02222 0.007550 0.006492 14.01 
0.02777 0.007650 0.006672 12.79 
0.03330 0.007750 0.006826 11.92 
0.05000 0.007900 0.007199 BAT 
0.08333 0.008350 0.007724 7.50 
0.13333 0.008650 0.008264 4e4T 
0.21667 0.009050 0.008879 1.89 
0.30000 0.009350 0.009323 0.29 
0.38333 0.010050 0.009674 3074 
063333 0.010500 02010430 0-66 
1.80000 0.011800 0.012122 2.73 
3¢30000 0.012900 0.013128 ~1.77 
3-80000 0.013050 0.013361 -2.38 
4.80000 02013500 0.013741 -1.79 
6203330 0.013900 0.014108 1.49 

10.08300 0.014600 0.014904 2209 
21464999 0.015900 0.016035 -0.85 
23264999 0.015900 0.016165 -1.67 
27216699 0.016250 0016371 -0.74 
28. 79999 0.016400 0.016459 -0.36 
30.04999 0.016400 0.016523 -0.75 
34229999 0.016650 0.016727 0246 
45263300 0.017200 0.917193 0.04 
48.01669 0.017300 0.017281 0.11 
50.01669 0.017400 0.017353 0.27 
51-86699 0.017500 0.017418 0047 
53284999 0.017550 0.017487 0.36 
58.20000 0.017600 0.017633 -0019 
69293300 0.018000 0.018003 -0.02 
71.83299 0.018000 0.018061 ~0234 
73259999 0.018050 0.018114 -0.36 
75493300 0.018100 0.018184 -0246 
77296700 0.018200 0.018243 -0.24 
61.66699 0.018250 0.018351 -0.55 
93043300 0.018550 0.018685 ~0.73 


C.5000C 
4.50000 
8.50000 
12.50000 
16.50000 
2C. 50000 
24.50000 
28.50000 
32250000 
36450000 
40.50000 
44.50000 
#8.50000 
52250000 
56-50000 
60.50000 
64.50000 
68.50000 
72.5000C 
76.50000 
80.5000C 
84-50000 
88.5000C 
$2.50000 
$6.-50006 


10. 
12. 
14. 
16. 
18. 
20. 
22. 
24e 
266 
28. 
30-6 
326 
34. 
36. 
38-6 
40. 
42. 
44. 
466 
48. 
50. 
526 
54. 
566 
58. 
60. 


1141.0 
1167.5 
1171.1 
1172.4 
1172.9 
1173.1 
1173.3 
1173.3 
1173.4 
1173.5 
1173.6 
1173.7 
1173.8 
1173.9 
1174.0 
1174.2 
1174.3 
1174.5 
1174.7 
1174.9 
1175.1 
1175.3 
1175.5 
1175.7 
1176.0 


1065.2 
1067.0 
1068.5 
1€69.9 
1071.1 
1072.3 
1073.4 
1074.3 
1075.3 
1076.1 
1077.0 
1C77.8 
1078.5 
1079.2 
1079.9 
1080.6 
1081.2 
1081.8 
1082.4 
1083.0 
1083.5 
1084.0 
1084.5 
1085.0 
1085.5 
1C86.0 
1086.4 
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[AL OF THE COMPOSITE 4NO TH 
O TIMES COVERING THE QURATI 
LEDASY IM ERED ERAS ON 
0.417 0.583 
0.403 0.597 
0.402 C.598 
0.401 0.599 
0.401 C.599 
0-401 0.599 
0.409 0.600 
0-400 0.600 
0.400 0.600 
0.400 0.600 
0-400 0.600 
0.400 0.600 
0.400 C.600 
0.400 0.600 
0-400 G.600 
0.400 0.600 
0.400 0.600 
0.400 0.600 
0.400 C.600 
0-409 C.6CO0 
0.400 0.600 
0.399 0.4601 
0.399 0.601 
0.399 0.6C)] 
0.399 C.6Cl 
CH 
ARE LISTEO 8ELOW: 
rEB*EY!M.  ERARLEBAGY HDA. 
0.460 0.540 
0.458 0.542 
0.457 0.543 
0.456 0.544 
0-455 0.545 
0.454 0.546 
0.453 0.547 
04453 0.547 
0.452 0.548 
0.452 0.548 
0.451 0.549 
0.451 C.549 
0.450 0.550 
0.450 0.550 
04449 0.551 
0.449 0.551 
0.449 0.551 
0.448 0.552 
04448 C2552 
0-448 0.552 
0.447 0.553 
04447 0.553 
0.447 0.553 
0.446 0.55% 
0-446 0.55% 
04446 0.554 
0.446 0.554 
0-445 0.555 
02445 C.555 
0.445 0.555 


U 


CORRESPONDING 
OF THF TEST 


+ 


626 
64 
66-6 
68. 
70. 
T2e 


4. 
Tbe 


78. 

80. 

82. 

84-6 

B66 

88. 

90. 

926 

94. 

966 

98. 
100. 
102. 
104. 
106. 
108. 
110. 
ll2. 
114. 
116. 
118. 
120. 
i22. 
124. 
126. 
128. 
130. 
132. 
134. 
136. 
138. 
140. 
142. 
144. 
146-6 
148. 
150. 
152-6 
154. 


156. 
158. 
160. 
162. 
164%. 
166. 
168. 
170. 
172. 
174. 
176. 
178. 
180. 
182. 
184. 
186. 
188. 
190. 
192. 
194. 
196. 
198. 
200. 


909.2 
908.8 
908.3 
907.9 
907.5 
907.0 


906.6 
906.3 


905.9 
905-5 
90561 
904.8 
904.4 
904.1 
903.7 
903.4 
903.61 
902.7 
902.4% 
902.1 
901.8 
901.5 
901.2 
900.9 
900.6 
900.3 
90061 
899.8 
899.5 
899.2 
899.0 
898.7 
898.5 
898.2 
898.0 
897.7 
897.5 
897.2 
897.0 
896.8 
896.5 
896.3 
896.1 
895.9 
895.6 
89544 
895.2 


855.0 
894.8 


894.6 
894.4 
894.2 
894.0 
893.8 
893.6 
893.4% 
893.2 
893.0 
892.8 
892.6 
892.4% 
892.3 
892.1 
891.9 
891.7 
891.5 
891.4% 
891.2 
891.0 
890.9 


1086.9 
1087.3 
1087.7 
1088.2 
1088.6 
1088.9 


1089.3 
1089.7 


1090.1 
1090.4 
1090.8 
1091.1 
1091.5 
1091.8 
10921 
1092.5 
1092.8 
1093.1 
1093.4 
1093.7 
1094.0 
1094.3 
1094.5 
1094.8 
1095.1 
1095.4 
1095.6 
1095.9 
1096.2 
1096.4 
1096.7 
1096.9 
1097.2 
1097.4 
1097.6 
1C97.9 
1098.1 
1098.3 
1098.6 
1098.8 
1099.0 
1099.2 
1099.4 
1099.7 
1099.9 
1100.1 
1100.3 


1100.5 
1100.7 
1100.9 
1101.1 
1101.3 
1101.5 
1101.6 
1101.8 
1102.0 
1102.2 
1102.4 
1102.6 
1102.7 
1102.9 
1103.1 
1103.3 
1103.4 
1103.6 
1103.8 
1104,0 
1104.1 
1104.3 
1104.4 
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0.445 
0.444 
0.444 
0.444 
0.444 
06444 


0-443 
0.443 


0.443 
0.443 
0.443 
0-442 
0.442 
0.442 
0.442 
0-442 
0.442 
0.441 
0-441 
0.441 
0.441 
0.441 
0.441 
0.441 
0.440 
0.440 
0.440 
0.440 
0.440 
0.440 
02440 
0.439 
0.439 
0.439 
0.439 
0.439 
0.439 
0.439 
0.439 
0.439 
0.438 
0.438 
0.438 
0.438 
0.438 
0.438 
0.438 


0.438 
0.438 
0.437 
0.437 
0.437 
0.437 
0.437 
0.437 
0.437 
0.437 
0.437 
0.437 
0.436 
0.436 
0.436 
0.436 
0.436 
0.436 
0.436 
0.436 
0.434 
0.436 
0.436 


0.555 
0.556 
0.556 
0.556 
0.556 
0.556 


0.557 
0.557 


0.557 
0.557 
0.557 
0.558 
0.558 
0.558 
0.558 
0.558 
0.558 
0.559 
0.559 
0.559 
0.559 
0.559 
0.559 
0.559 
0.560 
0.560 
0.560 
0.560 
0.560 
0.560 
0.560 
0.561 
0.561 
0.561 
C.561 
0.561 
0.561 
0.561 
0.561 
0.561 
0.562 
C.562 
0.562 
0.562 
0.562 
0.562 
0.562 


0.562 
0.562 


0.563 
9.563 
C.oes 
0.563 
0.563 
0.563 
0.563 
0.563 
0.563 
0.563 
0.564 
0.564 
C.564 
0.564 
0.564% 
0.564 
0.564 
0.564 
0.564 
C. 564 
C2564 


COMPOSITE at 1000 psi 
48.9% PE 


POLYETHYLENE at |000 psi 


a 
a 
O 
oO 
= 
a 
3 
LJ 
= 
LJ 
=! 
a 
a. 
oO 
© 
a. 
> 
eal 
© 
ou 


X-SCALE=4.Q0E-03 UNITS INCH. 
Y-SCALE=1.00E€+01 UNITS INCH. 


Ta et to rb Ie atte lates Ye LOU | ee) eG ete 
CURVE AND ACTURL DATA FOR COMPOSITE AB 
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%-SCALE=2.Q0E-O1 UNITS INCH. 
Y-SCALE=2.00E+01 UNITS INCH. 


HERING,BOX 21, LGAD FRACTIONS IN ERCH COMPONENT 
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APPENDIX D 
PREDICTION OF STRESS RELAXATION BEHAVIOR FROM EMPIRICAL CREEP EQUATION 
A stress relaxation test is performed by rapidly deforming a material 
to a predetermined strain then observing the decrease in stress as a func- 
tion of time while this strain is maintained. A series of such tests were 
run on the Instron testing machine for the plastics used in this study. 
The empirical equation developed for creep behavior was used to predict 


the relaxation curve. 


aya 
(D.1) a et Bott + eo (cl — ) 


Again the equation of state assumption was made for predicting relaxation 
behavior. Thus, for a given strain, and a series of times, the corres- 
ponding stress may be found by implicitly solving equation (D.1). This 
was done by a Fortran IV computer program called program RELAX. The input 
data consisted of the creep parameters previously determined and the ini- 
tial strain. A Newton-Raphson iteration was used to solve equation (D.1) 
for stress as a function of time, and these values were tabulated as out- 
put. <A plotting subroutine available at the NPGS Computer Facility (SUB- 
ROUTINE DRAW) was used to display the predicted relaxation curves. Fig- 
ure 32 illustrates a computed relaxation curve compared with the experi- 
mental relaxation curve during the initial part of the test. Figure 33 
shows the predicted curve and experimental data over the entire time frame 
of the experiment. These results are typical of all runs; the same general 
trend was observed for each plastic, even at low stress levels. The 
empirical creep equation predicts a greater decrease in stress than that 
observed experimentally. Apparently there is a structure factor influence 
creating irreversibilities that are unaccounted for in the creep equation. 


The equation of state assumption is not valid -- the stress level depends 
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on the strain history. 


An interesting observation can be made from the experimental data 
in Figure 32. During the first ten seconds of the test, the stress fell 
off rapidly in exponential fashion. This indicates the presence of mole- 
cular mechanisms with relaxetion times on the order of several seconds. 
The corresponding retardation cime in the creep equation is approximately 
one hour; the predicted curve does not show this sudden decrease at the 
start of the test. 't can be concluded, therefore, that the "instantaneous 


strain", parameter A, in the creep equation, includes a portion of the de- 


ii 
layed strain with very short time constants. The time scale of the creep 


test and the method of extrapolation to determine parameter A, prevents 


iL 
the resolution of this rapid but delayed elastic response. From Appendix 
A, Table 3, it is seen that the stress exponent N for the instantaneous 
strain is slightly greater than 1. If the delayed strain component com- 
bined in the term Ko™ couid be removed, this exponent would decrease. 
Although not a recessary consequence of the theory of rubber-like elasti- 
city [2], it appears that the exponent may approach a value of 1, so that 
under instantaneous loading the plastics behave as ideal Hookean solids. 

A modified equation was used with program RELAX, separating the para- 
meter A, into an instantaneous response and a rapid delayed elastic re- 


I 


sponse: 
(D. 2) E= Ko 
where 


(D.3) K a K , a Ka 


The stress coefficients Ky and K, were estimated from the experiment- 


al relaxation data and the time constant ae taken to be 10 seconds. 


IN 


a K oN ate 


The predicted relaxation curve then matched the data more closely at short 
times, but the same underestimate of stress level developed at long times. 


(See Figure 33). 
90 
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TABULATED OUTPUT FROM "RELAX" 


THE SPECIMEN AREA 180.04992 SQUARE INCHES 


THE CALCULATEO STRESS RELAXATION BEHAVIOR FOR NYLON IS LISTEN BELOW: 
TIME, HOURS TLME »SECONOS STRESS, PSI LOAO 
0.0 0.0 2000.0 99.8 
0.000278 le 1999.1 99.8 
0.000556 2. 1998.3 99.7 
0.000833 36 1997.5 99.7 
0.001111 4. 1996.7 99.7 
0.001389 5. 1996.0 99.6 
0.001667 Bhs 1995.2 99.6 
0.001944 7 1994.5 99.6 
0.002222 8. 1993.7 99.5 
0.002500 9. 1992.9 99.5 
0.002778 10. 1992.2 99.4 
0.003056 ll. 1991.5 99.4 
0.003333 12. 1990.7 99.4 
0.003611 13. 1990.0 99.3 
0.003889 14. 1989.3 99,3 
0.004167 15. 1988.5 99.3 
0.004444 16. 1987.8 99,2 
0.004722 17. 1987.1 99.2 
0.005000 18. 1986.4 99,2 
0.005278 19. 1985.7 9961 
0.005556 20. 1984.9 9961 
0.005833 216 1984.2 99.0 
0.006111 22. 1983.5 99.0 
0.006389 23. 1982.8 99.0 
0.006667 24. 1982.1 98.9 
0.006944 25. 1981.4 98.9 
0.007222 26. 1980.7 98,9 
0.007500 27. 1980.0 98.8 
0.007778 28. 1979.3 98.8 
0.008056 29. 1978.6 98.8 
0.008333 30. 1977.9 98.7 
0.008611 31. 1977.2 98.7 
0.008889 32. 1976.5 98.7 
0.009167 33. 1975.8 98.6 
0.009444 34, 1975.2 98.6 
0.009722 35. 1974.5 99.6 
0.010000 36. 1973.8 98.5 
0.010278 37. 1973.1 98.5 
0.010556 38. 1972.4 98.5 
0.010833 39. [e717 98.4 
0.011111 40. 1971.1 98.4 
0.011389 41. 1970.4 98.4 
0.011667 42. 1969.7 98.3 
0.011944 43. 1969.0 98.3 
0.012222 44. 1969, 98.3 
0.012500 45. 1967.7 98.2 
0.012778 466 1967.0 98.2 
0.013056 47. 1966.4 98.2 
0.013333 48. 1965.7 98.1 
0.013611 49. 1965.0 98.1 
0.013889 50. 1964.4 9A 
0.014167 51. 1963.7 98.0 
0.014444 52. 1963.1 98.0 
0.014722 53. 1962.4 98.0 
0.015000 54. 1961.7 97,9 
0.015278 55.6 1961.1 97.9 
0.015555 56. 1960.4 97.9 
0.015833 57. 1959.8 97.8 
0.016111 58. 1959.1 97,8 
0.016389 59. 1958.5 97.8 
0.016667 60. 1957.8 97.7 
0.019444 70. 1951.4 97.4 
0.022222 80. 1945.0 97.1 
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0.025000 
0.027778 
0.030555 
0.033333 
0.036111 
0.03AA89 
0.041667 
0.044444 
0.047222 


0.050000 
0.052778 


0.055555 
0.058333 
0.061111 
0.063889 
0.066667 
0.069444 
0.072222 
0.075000 
0.077778 
0.080555 
0.083333 
0.086111 
0.088889 
0.091666 
0.094444 
0.097222 
0.100000 
0.102777 
0.105555 
0.108333 
0.111111 
0.113888 
0.116666 
0.119444 
0.122222 
0.124999 
0.127777 
0.130555 
0.133333 
0.136110 
0-.13AAAR8A 
0.141666 
0.144444 
0.147221 
0.149999 
0.152777 
0.155555 
0.158332 


0.161110 
0.163848 


0.166666 
0.169443 
0.202777 
0.236110 
00269443 
0.302777 
0.336110 
0. 369443 
0.402777 
02436110 
0.469443 
0.502777 
0.536110 
0.569443 
0.602776 
0.636110 
0.665443 
0.702776 
0.736110 
0.769443 
0.802776 
0.836110 


90. 
100. 
110. 
120. 
130. 
140. 
150. 
160. 
170. 


LAO. 
190. 


200. 
210. 
220. 
230. 
240. 
.250. 
260. 
270. 
280. 
290. 
300. 
310. 
320. 
330. 
340. 
350-6 
360. 
370. 
380. 
390. 
400. 
410. 
420. 
430. 
440. 
450. 
460. 
470. 
4A0. 
4°0. 
500. 
510. 
520. 
530. 
540. 
550. 
560. 
570. 
580. 
590. 
600. 
610. 
730. 
ASO. 
970. 
1090. 
1210. 
1330. 
1450. 
1570. 
1690. 
1810. 
1930. 
2050. 
2170. 
2290. 
2410. 
2530. 
2650. 
2770. 
2890. 
3010. 


ao 


19348.A 
1932.7 
LaZ ha 
1970.7 


~£914.A 


1999.1 
1903.3 
1897.7 
1892.2 


1AR6.7 
1781.2 


1875.9 
1A7N,.4 
1865.3 
1840.2 
1855.1 
1850.0 
1845.0 
1840.1 
1835.2 
1830.3 
1825.8 
1820.A 
1A14.} 
1811.5 
1806.9 
1802.3 
1797,A 
1793.4 
1799.0 
1784.6 
1780.2 
1776.0 
L771. 7 
1767.5 
1763.2 
1759.7 
TIE It 
LSet 
1747.0 
1743.9 
1739.1 
1745.! 
1731.3 
1727.4 
1723.6 
1719.8 
1714. 
1712.3 
1708.4 
1705.0 
1791.2 
1697.7 
1656.7 
1619.1 
1584.4 
1552.9 
oa ee 
1495,7 
1469.9 
1445,A 
1423.1 
1401.7 
1381.5 
1362.4 
1344.3 
1327.1 
1310.7 
1295.0 
1280.1 
1265.9 
1252.2 
1239.2 


Q4,R 
96.5 
96.2 
95.9 
95.6 
95.3 
95.0 
94,7 
94.4% 


Q4.7 
93.9 


93.6 
93.4 
93.1 
92.9 
92.6 
92.3 
2a 
91.8 
91.6 
91.4 
91.1 
90.9 
90.7 
90.4 
90.2 
99.0 
AQ.7 
89.5 
89.3 
89.1 
88.9 
28.6 
BAG 
RAL? 
PAA 
A7,A 
R7.A 
R764 
R7.2 
RIN 
86,A~A 
R4K.6 
R6.% 
AB.? 
P60 
R5.,A 
R&,7 
aS.5 
A5S.3 
85.1 
or) 
R47 
A2.7 
80.8 
79.1 
Tee 
76.0 
74.7 
73.4 
Fee2 
71.0 
70.0 
69.0 
68.0 
67.1 
652 
65.4% 
64.6 
63.9 
63.2 
62-5 
41.9 


0 6869443 
0.902776 
0.936109 
0.969443 
1.002775 
1.036108 
1.536108 
2.036108 
22536108 
32036108 
32536108 
4.036108 
4.536108 
5.036108 
5.536108 
6.036108 
6.536106 
7.036108 
7.536108 
8.036108 
82536108 
9.036108 
9.536108 
10.036108 
10.536108 
11.036108 
11.536108 
12.036108 
12.536108 
13.036108 
13.536108 
14.036108 
14.536108 
15.036108 
15.536108 
16.036102 
16.536102 
17.036102 
17.536102 
18.036102 
18.536102 
19.036102 
19.536102 
20.036102 
20.536102 
21.036102 
21.536102 
22.036102 
22-536102 
23.036102 
232536102 
24.036102 
242536102 
25.036102 
25.536102 
26.036102 
26.536102 
27.036102 
27.536102 
28.036102 
282536102 
29.036102 
29.536102 
30.036102 


3130. 
3250. 
3370. 
3490. 
3610. 
3730. 
5530. 
7330. 
9130. 

10930. 
12730. 
14530. 
16330. 
18130. 
19930. 
21730. 
23530. 
25330. 
27130. 
28930. 
30730. 
32530. 
34330. 
36130. 
37930. 
39730. 
41530. 
43330. 
45130. 
46930. 
48730. 
50530. 
52330. 
54130. 
55930. 
57730. 
59530. 
61330. 
63130. 
64930. 
66730. 
68530. 
70330. 
72130. 
73930. 
75730. 
77530. 
79330. 
81130. 
82930. 
84730. 
86530. 
88330. 
90130. 
91930. 
93730. 
95530. 
97330. 
99130. 

100930. 

102730. 

104530. 

106330. 

108130. 
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1226.6 
1214.6 
1203.0 
1191.9 
1181.2 
1170.8 
1050.9 
972.3 
916.7 
875.3 
843.5 
A184 
798.2 
TAL.T 
7T6R.2 
756.8 
747.3 
739.3 
73264 
726.6 
721.5 
717.2 
713.8 
710.2 
T0764 
7046S 
702.7 
700.9 
699.2 
697.7 
696.4 
695.2 
694.2 
693.23 
692.5 
691.7 
691.1 
690.5 
689.9 
689.4 
689.0 
688.6 
688.2 
687.8 
687.5 
687.2 
686.9 
686.6 
686.4 
686.1 
685.9 
685.7 
685.5 
685.2 
685.0 
684.8 
684.7 
684.5 
684.3 
684.1 
683.9 
683.8 
683.6 
683.4 


61.2 
60.6 
60.9 
§9.5 
§9.0 
5A.4 
52.5 
48.5 
45.8 
43.7 
42.1 
40.9 
39.8 
39.0 
3A.3 
37.8 
37.3 
36.9 
326.6 
36.3 
36.0 
35.8 
35.6 
35.5 
35.3 
35.2 
35.1 

35.0 
34.9 
34,8 
34.8 
34.7 
34.7 
34.6 
34.6 
34.5 
34.5 
34.5 
34.4 
34.4% 
3404 
34.4 
34.4% 
34.3 

34.3 

34.3 
34.3 

34.3 

34.3 

34.2 

3402 

34.2 

34.2 

34.2 

34.2 

342 

34.2 

34.2 
34.2 
34el 
34.1 
34.1 
34.1 
34.1 


APPENDIX E 

PREDICTION OF STRESS-STRAIN CURVES USING THE EMPIRICAL CREEP EQUATION 

A simple method for predicting creep curves from stress-strain curves 
is outlined in Figure 34. As indicated in Figure 34(a), tensile tests are 
run at a series of constant strain rates. A cross-plot, Figure 34(b), is 
constructed for reciprocal strain rate vs. strain at a constant stress. 
By the graphical integration shown in Figure 34(b), times corresponding to 
certain strains may be found for a given stress level, and creep curves 
plotted as illustrated in Figure 34(c). This method implies the existence 
of an equation of state relationship between strain, stress, and strain 
rate, but does not explicitly state such an equation. A similar scheme 
was proposed by Marin and Pao [30], who obtained good agreement from test 
data on polymethylmethacrylate. 

In this study, the empirical creep equation was used to predict stress- 
strain curves by reversing the procedure outlined in Figure 34. From the 
empirical creep equation, 


_t® 
1 E€= KoN + Bo®t + PoM(i-2~ } 


The slope (strain rate) at any time for a given stress may be found by: 


-ct® 


@Q-| 
(E.2) ae = Bo" Per cat 


For a specified strain rate and a given stress, equation (E.2) may be 
solved implicitly for time, and this result substituted in equation (E.1) 
to obtain the corresponding strain. A Fortran IV computer program, called 
SIGEPS, was used to perform these calculations. Output from SIGEPS con- 


sisted of stress-strain data for a material for a given strain rate. 


of 


(a.) 


(b) 
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Figure 34 


Prediction of Creep Curves from Stress-Strain Data 
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Program SIGEPS and typical output is given on the following pages. Figure 
35 displays this output along with experimental tensile test data obtained 
at the same strain rate. The same characteristic was observed for most 
materials tested; the predicted tensile strain is greater than the observed 
strain for a given stress level. The deviation from the predicted curve 
becomes greater at higher stress levels. In the tensile test, stress 
changes with time, and internal structure changes probably develop. These 
changes are not accounted for in the creep equation used for the predic- 
tion. In addition, it should be noted that the tensile tests took place 
over a time span of about 10 minutes whereas the time duration of the creep 
test was 100 hours. Consequently the prediction of stress-strain curves is 
based on the creep curve at very short times. As indicated in the tabulat- 
ed output from program FITDATA (See Appendix B), it is at these short times 
where poorest fit of equation (E.1) usually occurred. 

Figure 36 illustrates an exception to the previous discussion. The 
predicted stress-strain curves for polycarbonate showed excellent agree- 
ment with tensile test data. It is possible that structure factor effects 
in polycarbonate are negligible at the stress levels encountered. Further- 
more, output from FITDATA for polycarbonate showed good agreement between 
the empirical creep equation and the test data at times under 10 minutes. 

It would be useful to further examine the feasibility of this method 
by utilizing an empirical creep equation specifically designed to give 


good data fits at short times. 
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Predicted Stress-Strain Curve for Nylon 
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Figure 36 


Strain , in/in x10? 
Predicted Stress-Strain Curve for Polycarbonate 
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TABULATED OUTPUT FROM "SIGEPS" 


POLYCARBONATE FOR VARIOUS STRAIN RATES IS LISTED BELOW 


THE STRESS=STRAIN OATA FOR 


IN/ IN/MIN 


STRAIN RATES, 


STRESS 


0.013040 


0.006087 


2000. 


0.007790 


2500. 


0.009535 


3000. 


0.011291 


3500. 


0.013212 


4000. 


0.015322 


4500. 


0.018104 


5000. 
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APPENDIX F 
THE BOX DISTRIBUTION FUNCTION FOR RETARDATION SPECTRA 

The use of a continuous distribution of retardation times to character- 
ize delayed elastic response was discussed in Section 6. Not only do 
theoretical considerations favor this approach, but mathematical solutions 
to viscoelastic problems may be considerably simplified, particularly if 
the functions representing the spectra can be handled by transform techni- 
ques. As an example, an analytic expression for creep behavior will be 
derived here, using a box distribution of retardation times as illustrated 
in Figure 37. 

An expression for creep strain corresponding to the four-element model 


of Figure 10, page 25 , may be written: 
-bt 
(F. 1) Ee= A, + Azt + Az (\- @ 5 


where the constant bb is the rate constant equal to the reciprocal of the 


characteristic relaxation time 7 , and coefficients A Ay» and An, de- 


lbs 
pend on the material and the applied stress. For a material with many dis- 


creet retardation times, the expression for creep strain given by equation 


(F.1) becomes: 


(F.2) e= A, +A,t + A, ma: (1-27) 
where the constant b; is the rate constant associated with the Le 
retardation time AR , and the factors Qi are weighting constants in- 
dicating the relative contribution of the im mechanism to the delayed 
strain. Since the factors QM are chosen so that 7,ai = \ » equation 


(F.2) may be rearranged as: 
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a(b) 





Figure 37 


Box Distribution of Rate Constants 


yas b 


Figure 38 


Wedge Distribution of Rate Constants 
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Ye Fe > 


3) €= Ay + At + Asl(i- Fae) 


If the delayed elastic response ( €qg ) has associated with it a continu- 


ous spectrum of retardation times, then the factors Aj may be represent- 


ed by a rate constant distribution function, q(b) , such that: 
b= 9 


(F.4) f a(b)db = | 
b 


=6 
For this continuous distribution of retardation times the creep equation 


(F.3) becomes: 


sy € = Ayvt Act + As(\- , q(b) 2° db) 


Equation (F.5) now contains the Laplace transformation of the rate con- 


Stant distribution function a(b) » and may be rewritten as: 


ws) € = A+ At + AZ(I- £ i ¢(b)} ) 


where 14 9(b) represents the Laplace transform of the function. 


The transform {iq (b)4 can be evaluated readily in the case of 


a box distribution of rate constants. Referring to the notation used 


for the box distribution illustrated in Figure 37, 


Lf 4b 


-— bt 
{4 2X db 
par 


oe [aoe Me 


7 Vy, gt b=c, 


il 





b=C, 


(F.7) Lt (by / ( ia = co ; 


In equation (F.7),constants Co and C, are the maximum and minimum rate 


constants, respectively, of the box distribution. These constants may 
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be evaluated for a given set of creep data as shown in Figure 19. 

In equations (F.1) and (F.6), the coefficient A. is intended to 
represent the delayed strain after elastic equilibrium is achieved, i.e., 
after infinite time. The delayed strain at time zero is zero. Therefore, 


the expression for {$4(b)} given by equation (F.7) must be normalized 


so that: 


@.s) —eligeiiaege , te=0 
9) ffa(by} 


This may be accomplished by expanding the exponential terms of equation 


i 
r- 
ck 

il 
8 


(F.7) in a MacLaurin series: 


fair} = tI! ct HYG. (gt at...) 


2. 2 
a Gah 


(F.10) 


; dim. ffa(byt = — Cc 


Equation (F.7) may now be normalized to get: 


| ie, cts 
eu L.{4lb} = Goat | ee | 


By substitution of the time limits ct =O and t =o into equation 
(F.11) it is apparent that the conditions of equations (F.8) and (F.9) 
are satisfied. Substitution of equation (F.il) into equation (F.6) gives 
the analytic expression for creep strain when the retardation spectrum is 


a box distribution of rate constants. 


LQ eels ~cot 
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Similar derivations may be performed whenever the retardation spectrum 

can be represented by a sum of functions whose transforms can be found. 
Unfortunately, the resulting expression is quite complex for all but the 
simplest distribution functions. For example, the wedge distribution func- 
tion, Figure 38, similar to the theoretical model of Alfrey discussed in 


Section 6, leads to the expressicn: 


eM" (kt-kt-i) +e" 


(F.13) €= A,+A,t + As| y Be he 


In the equations for creep strain developed in this Appendix, the 
delayed elastic response, Eg , represents a portion of the total strain. 


In equation (F.1), for example, 


-bt 
wis) €4 = A,(\-27 ) 


The value of the strain €g varies from zero at t= O to a value A, 
at t=, It is often convenient to normalize the delayed strain. The 


normalized form, called the delayed strain fraction, is defined by: 


(F.15) C= a io 


The delayed strain fraction varies from E=O at re ue to E€=\ at 
t = &© , and its value at any time t > O represents the progress 
toward elastic equilibrium. For equation (F.1), the delayed strain frac- 


tion is: 


—bt 
(F.16) aa he 


A better insight into the nature of the delayed elastic response is 
achieved by plotting delayed strain fraction versus the logarithm of time. 
Equation (F.16) is plotted in this manner in Figure 39. The inflection 


point of this sigmoidal shaped curve is reached at a time te = Yo , the 
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characteristic retardation time. The delayed strain fraction at t= VA, 

im E= \- Vo = 0.632. The slope of this curve changes very gradually 

in the range € = 0.3 to €=0.8. In this region, the shape of the curve 
can be approximated by a straight line drawn with a slope corresponding to 
AE | A (tnt) evaluated at €& = 0.5. This is shown by line segment A-A 
in Figure 39. 

For a box distribution of rate constants, the delayed strain fraction 
versus log time curve (Figure 40) is also sigmoidal in shape, with even 
less drastic slope changes in the region around €& =0.5. In order to 
accurately depict the delayed elastic response, the sigmoidal curve must 
fit the data closely at a time E, corresponding to € =0.5. Otherwise, 
large errors in predicted strain occur over most of the time range. Fig- 
ure 4l(a) illustrates a sigmoidal curve exhibiting a slope mismatch with a 
data plot. In Figure 41{b), the slopes of the data plot and the sigmoidal 
function are the same, but the responses do not occur at the same time. 

It is apparent from these figures that the error in predicted delayed 
strain fraction will be large if there is a slope and/or time mismatch at 
approximately ar 05. 

Plots of & versus aie were made for the creep data obtained in 
this investigation. For different viscoelastic materials tested, these 
plots exhibited remarkable similarities. The delayed strain fraction 
reached a value € = 0.5 atatime t ~1 hour, and the slope of the 
curves at E€ = 0.5 was AE Jol (Int) == 0.15. The same values were observed 
for PP, PC, NY, and PE and were independent of stress. Therefore, the 
same function may be used to characterize the delayed elastic response 
of each of these materials. A question arises as to whether or not a func- 
tion such as the box distribution leads to a suitable representation. To 
text the applicability of the box distribution, it is necessary to investi- 


gate the characteristics of the resulting sigmoidal curve for the limiting 
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O Int 


Figure 39 


Delayed Strain Fraction vs. Ln Time for a Single Retardation Time 





Int 


Figure 40 


Sigmoidal Curve for a Box Distribution 
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Unsuitable Representations of Delayed Elastic Response 
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cases of the box distribution shown in Figure 42. This analysis is given 
below: 
If the limits of the box distribution are extended to cover many de- 


cades of time, as shown schematically in Figure 42(a), 





a ey 
ee ~ Cae 
yes Eg FOr a K® 
The expression for the delayed strain fraction (from equation (F.12)) is 
given by: 
GHG -Czt 
g= \- | ££ — 
(Ex) Gey = Cet 
The terms gat and C, may be smerlected so that; 
wall 
zo ! 
F.18 = a 
( ) a | Cat 


The constant Cc, is very small so that, 


=e = | 


The delayed strain fraction for the extended box distribution may then be 


approximated by: 


(aie) eee rants Ce 7 O 


This curve is illustrated in Figure 43(a). The approximation is obvious- 


ly in error for times t <a Vo, » but the constant C, is large and the 


2 
equation (F.19) is valid in the time region of interest. 

For the case shown in Figure 42(b) where the limits converge and the 
box becomes narrow, the equation for delayed strain may be approximated 


as follows: Denoting the bracketed term of equation (F.17) by A >» and 


expanding the exponentials, 


al 7 2 3 | 2 P 3 


\ 


eS (C,-¢,)-t 
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(i. 20) f — (Gi=G) a aac) t ay Care) 3 
UIX< 1- Gs) pee) ay O) 


ee ee Oe Ge ansiOm 1G Kenresented by: 
fv" A 
Ce. 5 aE ad 


- 
~ o é 
‘ \ Lae f sa) 
em oe omer terme cace of converzing limits ( CC, -- = Gp 


2 st 


may be dectemined by eonsicering C, and C, to differ by a small quantity € , 
ns i 


8 


aC —s o —_ ia ot . 
a) ad Ge * > ea = 
then 
Ww - We 
CG Gea ) 


Using a binovial expansion and truncating after the second ter, 





A} 2 _\/ Sa "3 a | 4 
a eee) = Oe eC, , for €C&e, 
MA Aa tie Ane - 
then, - - ae Z en mee es Pa eiG ov \ 





mill," 6.) m\ (0. 4+0-¢,) ale (tr-sj! 
therefore, 
“ WN SA. a 
eon 7 es C, C 


(i) ee aes 


C2775, ml (tf, -0,) = Gein 


th 


Suvstituting etuaticn (F.21) for the m tema ef ecuation (F.20) and 


evaluatin; the torms for #7L= 1,2,3..... .the expansion becomes: 


if mM~! 
ec cet? Co M~} 
t 
(F e 2 e + = } aa ee -j a 42 t £ 6 ee ee 
m-1) 


<a 


| mr 


SOU 


Equation (F.22) reoresents the expansion of 92 in a MacLaurin 


series, so the desired approximation has been found: 


sabe 


ary 
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Equation (F.23) would result from a "single spike" distribution, corres- 
ponding to a single characteristic retardation time. The analysis veri- 
fies the expected result that equation (F.17) should reduce to equation 
(F.23) as the width of the box distribution becomes vanishingly small. 
This approximation is sketched in Figure 43(b). 

In order to be used for an accurate representation of the data in this 
investigation, the slope of the sigmoidal plot of the box distribution must 


take on the value: 
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For the case of narrowing box width, (Cc, —> C,)> given by equation 
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Figure 42 


Limiting Cases of the Box Distribution 
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Figure 43 


Approximations of the Response Function for the 
Limiting Cases of the Box Distribution 
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It may be demonstrated by numerical methods that equation (F.26) represents 
the minimum value of the slope AE | Unt) for any choice of the 
constants C, and C,s provided C, ? C,. The requirements of equation 
(F.24) can not be met, since the observed slope of 0.15 does not fall in 
the range of 0.35 to 0.5. Therefore the box distribution is an unsatis- 
factory choice for the retardation spectrum of the materials tested in 

this experiment. 

A similar analysis of the wedge distribution (Figure 38) shows a 
slope range 0.35 < (dé fd (Mt) |e =0,5) < 0.4 for all choices of KR and K,» 
indicating that this rate constant distribution is also unsatisfactory. 
Since the limiting case of the box distribution given by equation (F.23) 
corresponds to the response of a single Voigt Element (Figure 5), the four 
parameter creep model (Equation (F.1)) gives a poor representation of de- 
layed elastic strain. 

As discussed in Section 6, good correlations can be obtained using 
a model with several discreet retardation times. This "multi-spike distri- 
bution" produces a plot of delayed strain fraction versus log time with 
inflections at each retardation time. These retardation times may be chosen 
so that the response function approximates the creep data quite closely. 
Figure 44 illustrates the response function resulting from a three-spike 
distribution superimposed on a typical plot of creep data. The response 
function from a box distribution is also shown. The lower limit, C,> of 


the box is equal to the slow rate constant b, of the three-spike distri- 


bution, and the upper limit C, equals the fast rate constant bs 
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These values for the box distribution function were chosen to illustrate 
the dominant role of the fast rate constants in governing the delayed 
elastic response. The constants C, and C,, could be changed to shift the 
sigmoidal curve toward longer times, but the slope cannot be adjusted to 
produce a good fit with the data. 

Monomodal distribution functions characterized by a single maximum 
can always be approximated by a sum of box-type distributions as shown in 
Figure 45. The limits of each succeeding box fall within the range of the 
first box. Since the slope of the delayed strain fraction versus log 
time plot resulting from such a distribution is restricted in accordance 
with equations (F.25) and (F.26), it is apparent that the retardation 
spectrum of the plastics tested is not monomodal. 

Bimodal (or polymodal) spectra produce a series of inflections in the 
plot of €& vs. byt. » which results in a lower overall slope, i.e., the 
delayed elastic response occurs over many decades of time. The three- 
spike distribution of Figure 44 is an example of a polymodal retardation 
spectrum. 

The delayed elastic response of the materials tested covered about 
five decades of time. This response can be characterized by a polymodal 
retardation spectrum such as illustrated in Figure 46, where the peaks 
imply the presence of several different molecular mechanisms whose time 
constants differ by orders of magnitude. This model is in accord with 
Alfrey's concept (See Section 6) of local bond straightening and long- 


range uncoiling of the polymer chains. 
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Figure 49 


Monomodal Retardation Spectrum Approximated By Box Distributions 
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Figure 46 


A Polymodal Distribution Function 
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APPENDIX G 
EQUIPMENT 

The creep machines used in this investigation apply constant dead- 
weight loading through calibrated lever arms. (Figure 47). Eleven identi- 
cal units are installed in one test stand. Photographs of the laboratory 
and apparatus are shown in Figures 48 through 50. Specimens were gripped 
by steel pins placed through the 5/8 inch diameter hole in each end, and 
uniaxial loading was maintained through a series of universal joints 
above and below the grips. Turnbuckles and an adjustable counterbalance 
weight permitted final balancing of the lever arm prior to load applica- 
tion. 

Studies of viscoelastic behavior of materials require large amounts 
of data obtained over long time periods of testing. Consequently, facili- 
ties for automatic recording of creep and relaxation data are highly de- 
sirable. Design of such instrumentation was started in conjunction with 
this project and is presently being installed. Special clip-on strain 
gage fixtures were fabricated to permit measurement of the large strains 
encountered in mechanical testing of viscoelastic materials. These clip 
gauges, shown in Figures 51 and 52 consist of phosphor bronze strips with 
BLH SR-4 Type A-5-1 strain gauges mounted on each side. Tests on phosphor 
bronze indicate that this material has remarkable linearity over the en- 
tire elastic range with no observable hysteresis effects. The strips are 
connected to aluminum legs which make a knife-edge contact with the speci- 
men. This design has two useful features: 

(1) The legs are detachable so that different leg lengths may 
be used depending on the strain sensitivity desired. 

(2) The mode of attachment to the specimen is such that this 
electromechanical device may be used for specimens of various 


widths and thicknesses. 
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The four strain gauges are connected in a Wheatstone bridge arrange- 
ment shown in Figure 53. A potentiometer in the circuit permits a null 
adjustment to be made. This bridge circuitry has several advantages: [26] 

(1) Automatic moment compensation (in the event there are devia- 
tions from axial loading. 

(2) Signal strength increased by a factor of four. (over that 
of a single gauge mounted on the specimen) 

Terminal boards have been set up on the eleven-unit creep machine to 
receive the leads from the extensometer circuit for each test specimen. 
Bridge circuits, each with its own balancing potentiometer, are placed 
in parallel across a D.C. power supply, with leads from each bridge enter- 
ing master control panels. From these panels, the output signal from 
each extensometer may be fed into a multi-channel recorder, permitting 
testing and graphical data output for up to eleven specimens simultaneous- 
ly. A stress analysis of the clip gauges, using linear elastic theory, 
is given in Figures 54 and 55. 

It is desirable to have the output of electromechanical extensometers 
converted directly into digitalized form because of the usefulness of the 
digital computer in data analysis. Preliminary design of such instru- 


mentation has begun at NPGS, Monterey. 
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Lever Arm Loading Assembly 
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Specimen Grips and Clamps 
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Figure 53 


Bridge Circuitry for Extensometer 
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Dimensioned Sketch for Clip Gauge Stress Analysis 
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CLIP GAUGE STRESS ANALYSIS 


o Assume no bending in aluminum legs 
« To find E in gauges for a given specimen extension, 0 


Af 
po = 5 : a= ain~'Ge 
from the flexure formula, €= 4/0 


p= frzam(®) yas, 


es é ime (ti. 


for small angles ©, im' (2. = S 
e= *¥eLh 

for d=.0i0, h=2", L=? 
E= .00125 § 


for a specimen extension of § =}" 
E = 1250 Kinches/inch 
or about one-half the maximum strain 
capability of the BLH type a-5-l gages 


maximum stress in the bronze strip: 
_e €=.00its 
O=EEe, E= 1G x 10 psi 
O = LOOOO PSs\ 


or only a small fraction of the yield 
strength of phosphor bronze (Oy = 112,000 psi) 


at an extension of $= {" , error in the strain 
due to the small angle assumption is: 


Son = Vor =.25 , ain’ (3) =.2474 


025 —-.1474 2 iy 
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¢ Error = 75 


Figure 59 


Theoretical Stress Analysis of Clip Gauges 
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